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Abstract 

The response of particle density to a dilation of a periodic potential in an insulator, 
with or without a fixed background potential or a magnetic field, is shown to be quantized. 
A similar phenomenon occurs in a quantum Hall system, where the derivative of the 
electron density with respect to the magnetic field is quantized in units of e/h, the inverse 
flux quantum. A number of other interesting results may be unified under the notion of 
quantized density response, including the quantization of adiabatic particle transport, the 
gap labeling theorem and its generalization in the presence of a magnetic field. 
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1. Introduction 

It is well known that for a two dimensional gas of electrons in a magnetic field B, 
the particle density p satisfies p/(eB/h) = n for n filled Landau levels below the Fermi 
energy, where h/e is the flux quantum. In the presence of a weak periodic potential, each 
Landau level is split into p subbands, if the magnetic flux in each unit cell is | times the 
flux quantum, where q and p are mutually prime integers. When the Fermi energy lies in a 
gap between two such subbands, the ratio p/(eB/h) is no longer an integer, because each 
filled subband contains a particle density p times smaller than that for a full Landau level. 
What remains to be quantized is instead the derivative 



which is also the integer for the quantized Hall conductance. [1-4] Indeed, according to the 
Streda formula, [1] the Hall conductance may be generally written as 



whenever the Fermi energy p lies in an energy gap (or even a mobility gap). Under the 
same condition, the Hall conductance and therefore the density response are quantized. 
Note that disorder is certainly allowed under this condition. 

In this paper, we would like to show that the notion of quantized density response 
may be applied to other type of insulators. Consider now an electron gas in a periodic 
potential. The average particle density satisfies pv = n for n filled Bloch bands, where v is 
the unit cell volume. A natural question is whether one can generalize such a result to the 
case with an additional background potential and/or a magnetic field. The background 
potential may be periodic with a different periodicity or even be disordered. The magnetic 
field does not have to be uniform. In analogy with the quantum Hall system discussed 
above, we would like to expect that the density response to a dilation of a periodic potential 
should be quantized when the Fermi energy lies in a gap (or even a mobility gap), although 
the density itself may not be quantized. We show in this paper that this is indeed the 
case for a number of interesting situations. We would like to emphasize that the notion of 
quantized density response may provide a unifying framework in which several important 
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subjects, such as the quantum Hall effect, quantized adiabatic particle transport, and gap 
labeling theorem for quasi-periodic systems, are found as applications. 

The organization of our paper is as follows. Our initial study will be limited to 
the relatively simple case of one dimensional systems (Section 2), but the result will be 
generalized to two dimensions afterwards (Section 3). The presence of particle-particle 
interactions makes the problem complicated, but does not alter our conclusion at least for 
one dimensional systems (Section 4). Density response to a magnetic field will then be 
studied in Section 5. Some useful applications of our results will be discussed in Section 6. 

2. One Dimensional Case 

Consider a system of non-interacting electrons in the following potential at zero tem- 
perature: 

V(x) = V 1 (ax-0 + V 2 (x), (3) 

where V± is a periodic function of £ of period 1, so that the periodicity in x is 1/a. We 
will assume that the Fermi energy lies in an energy gap for arbitrary £ and for a small 
neighborhood of a. We would like to show that is quantized as integers, where p is the 
density averaged over position in the bulk of the system and over £ in a period (unity). 
The average over £ will be convenient for the mathematical manipulations, but is not really 
necessary if the potential V 2 is random, in which case there is no prefered position for V\. 
The same is also true when V 2 is periodic but incommensurate with V\. 
The particle density at x$ may be written as 

P (x ) = f^Tr{Gp(xo)}, (4) 

where G = (z — H)~ x is the resolvent of the single particle Hamiltonian, and p{xq) = 
S(x — xq) is the density operator. The integral goes along a contour in the complex energy 
plane enclosing the spectrum of filled states. The contour will be chosen everywhere away 
from the real axis where the spectrum lies, except at the Fermi energy and a point below 
the whole spectrum. We assume that there are no extended states in the bulk at the Fermi 
energy. It then follows that the Green function in the coordinate representation G(x,x') 
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is essentially local, i,e., it decays exponentially at large \x — x'\.[5] Our later manipulations 
will depend critically on this property. 

Using the differential formula dG = —GdG~ 1 G, the density response can be written 

as 



da J 2ni [ <9£ 

We have expressed the a-derivative of V\ in terms of its ^-derivative. Next, we note that 
xG may be written as x]G + Gx = GihvG + Gx, where v is the velocity operator. 

Using further the fact that dVl ^~^ = , we may write (5) as 

dp(xp ) 
da 

When we average over £ in a unit cell, the second term on the right hand side of (6) 
becomes zero because of the periodicity of G in £. In the following, we will no longer 
consider this term. 

The locality of the Green function allows us to do a number of things. First, the 
unboundedness of x should not introduce any difficulty in the above manipulations, if we 
measure the density only in the bulk of the system. Secondly, the locality of G makes sure 
that the expression (6) is insensitive to changes in the boundary condition. In particular, 
we may impose the quasi-periodic boundary condition ip(x + L) = e ie ip(x), and make an 
average over 9. Under such a boundary condition, the average of p(xo) over xq in the bulk 
may be achieved by integrating over xq in the whole space and dividing it by L. Therefore, 
with p also averaged over £, we have 

dp 1 f 2lT d6 f 1 r dz ( 



. — d£(b Tr ^ G — -\r- — ^GihvG } . (7) 

da L Jq 2tt Jo J liri \ d£ J 1 ' 

To show explicitly the quantization of the density response, we make a gauge transfor- 
mation ip(x) — > 4>(x) = tfj(x)e~ l0x / L , then G — > G, and v — > — -| 9< ^ e , and the expression 
(7) may be written as 

-27T jo r r a-, { Pkn-t a/^-i 



f 2n de f „ l dz ^ f ^dG- 1 - dG- 1 ~ 



(8) 



where the second line is derived from the first by using the fact that GG = — ^ and by 
integrating over z by parts. According to Ref.[6], the above integral takes integer values 
and represents a topological invariant for mappings from T 3 to nonsingular operators. 

Finally, we would like to remark that the condition that the Green function falls off 
exponentially at large distances is sufficient but not necessary. A sufficiently fast power 
law fall off (e.g. \x — x'\~ 3 in ID) should still be alright for the quantization of density 
response in the thermao dynamic limit. Such a remark also applies to situations to be 
discussed in the subsequent sections. 
3. Two Dimensional Case 

The generalization of our result to higher dimensions involves surprisingly a lot more 
algebra. We have so far only achieved this for two dimenions. The potential V\ in (3) now 
has the form Vi(a\Xi — £1, a 2 x 2 — £2), which is periodic in £1 and £2 with unit periodicity. 
We would like to show that 

da\da 2 ^ 
is quantized as integers, where p is the density averaged over the position in the bulk of 
the system and over £1 and £2- Again, the average over the £'s is not really necessary if 
V 2 is random or incommensurate with V\. We will also assume that there are no extended 
states in the bulk at the Fermi energy, but states at the boundary, extended or not, are 
allowed. The system may also be in a magnetic field, which may or may not be uniform. 
The locality of the Green function will again be the key to our problem. 

The general method is quite parallel to the one dimensional case. Equation (5) now 
becomes 

d 2 p(r ) r dz , u 97i dVi 

+G W, X2G W, XlG (10) 

The term including the double derivative of V\ may be rewritten as 

-G^G Xl A~ GxiX ^G^G, (11) 



where we have ignored a term of total derivative in £1, because such a term will become 
zero when we average over £1 in (9). The quantity in the square bracket in (10) may now 
be written as 

G— xi, G] -wrx 2 G + Gx 2 — - [G, xi j—G 

r dV x dV 1 dV x dV 1 { } 

=G——GihviG——x 2 G - Gx 2 ——Gihv 1 G——G, 
d£i d£ 2 d£, 2 d^i 

where we have used the fact that [xi, G] = x{\G = GiTiv\G. Next, we may replace 

x 2 G in the first term and Gx 2 in the second term of the last expression by ±[x 2 ,G] = 
±Gihv 2 G respectively, because the error made in (10) is proportional to a total derivative 
in £i and 

Tr{p(v )[G^Gv 1 G^Gx 2 - X2 G^Gv 2 G^G}} 

8 8C 8 8C 

= lV Wro ), 2[ _( Gci _)-_( Gci _)] } , 

where we have used the cyclic property of the trace and the commutativity of x 2 and /3(r ). 
As a result of these considerations, we may now rewrite the quantity in (10) averaged over 
£i and £2 as 

\ \ d^dfr <b -^Tr I p(r )[G—±Gih Vl G—±Gihv 2 G + Gihv 2 G—±Gih Vl G—±G] 
J J J 2m { <9£i d£ 2 d£ 2 d£i . 

(14) 

Finally, as in the one dimensional case, we introduce the quasiperiodic boundary conditions 
in both x and y directions, and make average over the phase angles 6>i and 9 2 , then the 
quantity in (9) becomes 



At At II d h^l I idz rr ^dG- 1 -dG- 1 -dG- 1 -dG- 1 ~ 

-dG^-dG- 1 -dG- 1 -dG- 1 ~, 

+G ^e 2 - G ^ G ^eT G ^ G ^ 



(15) 



where we have made a gauge transformation so that the phase angle dependence in the 
boundary conditions goes explicitly into the Hamiltonian. 

In the appendix we show that (15) can be further written as 

Mid^d^hdfhdz s-dG^-dG^-dG^-dG^-dG- 1 
iso^ Tr{G ^r G ^eT G ^ G ^e 2 - G ^- (16) 

+119 antisymmetrizing terms}. 



According to Ref.[7], such an expression takes integer values and represents the winding 
number for the mapping from T 5 to nonsingular operators G. Therefore, the density 
response (9) is quantized as integers. 

4. The Case With Many Body Interactions 

In this section we try to generalize the result of quantization of density response to the 
case with many body interactions. This can be readily done by invoking the relationship 
with quantized particle transport (see the next section), which has been established in the 
presence of many-body interactions. As we show below, such a generalization can also be 
done directly within the present formalism in one dimension. 

Eq.(5) remains valid if we regard G as the Green function for the many body Hamil- 
tonian, regard p(xo) as J2f=i^( x j ~ x o), and replace dVl ^~0 x D y J2f=i 9Vi (°^j — Xj, 
where j labels the particles. The contour integral is now surrounding the ground state 
energy of the many body system. Eq.(6) will also remain valid if we make similar changes 
on the right hand side; explicitly, the second term becomes 

/^E^{G^t^w}, (17) 

where the trace is over the antisymmetric many-body states only. In the coordinate rep- 
resentation, G(xi, X2, • • ■ , xn; x[, x' 2 , • ■ • , x' N ), the Green function is antisymetric under 
permutation of each set of position variables. We may thus replace Xj in the above equa- 
tion by ^2,Xk/N under the trace operation. Then Eq.(17) becomes a total derivative in £, 
and can therefore be ignored. Similarly, Eq.(7) now becomes 

dp(x ) N [ 27T d6 f I dz ^ / dV^aXj - o 



da L J 2n 



J d( f^yj^o GaVjG }. (18) 

Again, because of the antisymmetry of the Green functions in particle labels, we may 
replace Vj by J2 k Vk/N. The final result of (8) can then be rederived, showing the quanti- 
zation of density response in the presence of many body interactions. 

Unfortunately, similar arguments in the two dimensional case show that the density 
response is given by the expression (15) or (16) divided by N, if we regard G as the many- 
body Green function. This says that density response is only quantized in integer units of 
l/N. 



5. Density Response to a Magnetic Field 



Having studied the density response to the dilation of periodic potentials, we now 
consider the problem of density response to a change of a magnetic field. Associated with 
the magnetic field is a length, the so called magnetic length, whose square gives 

the area occupied by one flux quantum. Our problem can therefore also be stated as how 
the electron density is controled by the scale of the magnetic length rather than by the 
unit cell length scale of a periodic potential. As we pointed out before, this is a problem 
closely related to the quantum Hall effect, but our treatment will be shaped in a way more 
parallel to the method used in the previous sections than those used by other authors. [8] 

Consider a two dimensional electron gas in a uniform magnetic field B in the perpen- 
dicular direction. The single electron Hamiltonian may be written as 

# =^(p + eA) 2 + V(r) (19) 

where A is the vector potential and may be taken as (A\, Bx + A 2 ) with A\ and A 2 being 
constants independent of B. We have no special requirement on the scalor potential V(r). 
The particle density is still given by a formula like in (4) except that p{xq) should be 
replaced by the corresponding 2D quantity p(ro). Since the magnetic field dependence 
comes only from the Green function, the differential density response to the magnetic field 
can be calculated using the formula dG = —GdG~ l G as 



-/!H G !!«} (20) 



Next, we note as before that xG may be written as G[G~ X , x]G + Gx = ^G-j^-G + Gx. 
Equation (20) can then be written as 

d P (r ) h rdz ( dH dH \ 

-dB- = —eJ 27 Tr yd~A 2 d~A\ J (21) 

where we have dropped a term proportional to the ^-integral of Tr f^G-j^Gxp(r ) j. Such 
a term is proportional to x$ times the derivative of p(r ) with respect to A 2 , and is zero 
because of gauge invariance. 



As before, the z-integral goes along a contour in the complex energy plane enclosing 
the spectrum of filled states. The contour will be chosen everywhere away from the real 
axis (where the spectrum lies), except at the Fermi energy and a point below the whole 
spectrum. We assume that there are no extended states in the bulk at the Fermi energy. It 
then follows that the Green function in the coordinate representation G(r, r') is essentially 
local, i,e., it decays exponentially at large |r — r'|.[5] We can thus impose periodic boundary 
conditions to the wave functions, making only exponentially small errors in the density 
response in the bulk. We can then average over tq on the boundaryless torus, yielding 

M-o) ft I^lo^a^c) (22) 



dB eLiL 2 / 2tt [ dA 2 dA 1 

Also upto exponentially small errors, we can make a further average over A\ and A 2 in 
the intervals (0, ^-) and (0, -^) corresponding to unit flux quanta through the holes of 
the torus, with the result 



which is essentially the same as Eq.(l) of Ref.[6], because the above expression is antisym- 
metric upon interchanging A\ and A 2 , and because the ^-derivative of G~ x is unity. The 
density response to the magnetic field is therefore quantized in units of e/h. 

As has been shown in the context of the quantum Hall effect, the above result allows a 
number of generalizations. First, the system does not have to be strictly two dimensional. 
Also, one can add a modulation to the magnetic field without changing the quantization 
of the density response to the uniform part of the field. Further, one can reformulate 
the above derivations using many-body quantities (cf. Section 4), allowing interactions 
between electrons. Finally, when the many-body state does not go back to itself when a 
flux quantum is inserted through a hole of the torus, which occurs in a fractional quantum 
Hall regime, one obtains fractionally quantized density response to the field. 
6. Applications 

In one dimension, quantized density response to a dilation of a periodic potential is 
the same as the quantized adiabatic particle transport studied by Thouless and others. [8,5] 
When a is increased by 5a, the potential Vi shrinks a distance of 4^ at x = ±L/2. If 



the particle transport is t when the potential is shifted by a period (1/a), there will be 
tadaL/(2a) = tSaL/2 particles transported inward through the point x = L/2, and there 
will be the same amount of particles transported inward through the point at x = —L/2. 
The density increase is therefore tSa. In fact, the expression (7) is identical to that for the 
particle transport. However, the concept of density response can be readily generalized to 
more than one dimension. 

The quantized density response as a topological invariant has an important physical 
content: the system may be classified into topologically distinct situations, two situations 
being topologically the same if one can be continuously transformed into the other without 
closing the gap at the Fermi energy. Two situations should have the same integer of density 
response if they are topologically identical. This follows from the fact that the expression 
for the density response is a continuous functional of the Hamiltonian, and from the fact 
that the value of the density response is discrete. 

One special class is for V± = const, or weak, where we have to assume that V2 is such 
that there is always a gap at the Fermi energy. A V\ potential is considered weak as long 
as it does not close the energy gap at the Fermi energy as it is turned on from zero. The 
integer of density response for this class is zero. 

Another special class is for V2 = const, or weak. For this class we only need to consider 
the case of V2 = 0. The density response for n filled bands is simply n, because there is 
one electron per unit cell for each filled band, discounting the spin degeneracy. 

More interesting situations occur when the potential is a superposition of periodic 
ones. Consider the following example in one dimension: 



where each Vj has period 1 in £j. We assume that the gap at the Fermi energy does not 
close for a neighborhood of ctj, and that tj is the integer of density response to the jth 
potential, then 



3 

Under a uniform scaling of x — > (1 + Sa)x, we have Sp = pda. On the other hand, such 
a scaling is equivalent to setting Sctj = ctjSa, so that dp = da ^ ■ tjctj, according to (25). 




(24) 
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(25) 
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Comparison of the two results yields 

P = ^2 t j a 3- ( 26 ) 

3 

A similar result has been obtained by Thouless using quantized particle transport. [9] In 
higher dimensions, the present theory allows us to generalize the above result to 

P = Y,tjhj, (27) 

3 

where Vj is the unit cell volume of the jth potential. [10] 

In two dimensions and in the presence of a magnetic field B, one has to add to the 
right hand side of (27) an extra term aeSB/h, where a is the integer for the quantum Hall 
effect. Under a uniform scaling of da, we have to also set SB = SaB, then we have 

p = ^2tj/vj +aeB/h. (28) 

3 

Such a result was first derived by Wannier in the study of the Harper's equation[ll], and 
has been proved in general by Dana et al using magnetic translation group. [4] A three 
dimensional generalization of the above relation has been studied in Ref.[12]. 

There are two remaining problems to be solved. The first is to generalize the result of 
density response quantization to three (and higher) dimensions. It may be possible to do 
so within the present formalism, but the algebra involved will be quite heavy. The second 
is to show (or disprove) the quantization in two or more dimensions with particle particle 
interactions present. It is possible that fractionally quantized density response may occur 
for the interacting case. [13] 
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Appendix 



To simplify notations, we represent the quantity in the square brackets in (15) sym- 
bolically as 

The original problem is symmetric under the exchange of (1 2), (Al) is equivalent to 
half of 

£l#l£2#2 + #2£2#l£l 



+&02£i0i + 0i£i02&. 
By integrating z by parts, we can rewrite (Al) as 

-0i602fi-&02£i0i-02fi0i6 

— £2#l£l#2 — ^1^1^2^2 — 

Therefore, (A2) is equal to 

- 0i&0 2 £i - #2 W2 
It then follows that (Al) may be written as | of (A2) plus (A4), i.e. 

+02^1 - o 2 ZiOi& 

+^2^2^1^1 — £2#l£l#2 
+^1^1^2^2 — #l£202£l- 

The first line of (A5) may be written as 

~ dG' 1 dG dG' 1 dG _ ^dG' 1 dG dG' 1 dG 

dfi W x d£ x W 2 ~ dl[i d9 2 d^i df x 

dG dG- 1 - dG- 1 dG dG dG~ x ^dG~ x dG 
-G——— + — — — — G- 



dQ x dix dix d9 2 d6 2 dfi d^i ddx 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



where we have integrated 6>i by parts in the first term, and 9 2 in the second term. Written 
in abstract notation, (A6) is 6*2^1^2^1 _ #i£i£2#2- Similar manipulations can be made on 

12 



the remaining lines of (A5), showing that (A5) is equal to 

(A7) 

+£l#2#l£2 — £2#2#l£l- 

Next, by integrating z by parts, we can show that (A7) is further equal to 

£l£2#2#l — £l£2#l#2 
+£2£l#l#2 — ^2^l6*2^1 

048) 

+6*1^2^2^1 — ^1^2^1^2 
+#2#l£l£2 — ^2^1^2^1- 

In summary, we may write (Al) as ^ of (A5)+(A7)+(A8), namely 

-dG- 1 -8G- 1 -8G- 1 -8G- 1 ~ nn . 
G G G G G + 23 antisymetrizing terms, 
<9£i dd\ <9£ 2 o6 2 

where 'antisymmetrizing' acts on the four variables £i,0i,£2> an d #2- Then, we append 
each term a factor of d ^ z = 1, and use the cyclic property of the trace, yielding the 
expression (16) in the text. 
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